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Graphene nanoribbons are the flimsiest material systems in the world, and they get readily dis- 
torted. Distortion by twisting, for one, is important because it couples to ribbon's electronic prop- 
erties. In this Letter, using simulations with density-functional tight-binding and revised periodic 
boundary conditions, I show that twisting appears almost equivalent to stretching; electronic struc- 
tures in a given nanoribbon either upon twisting or upon certain stretching are quantitatively similar. 
This simple equivalence will provide a valuable guideline for interpreting and designing experiments 
with these flimsy ribbons. 
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The fiimsyness of graphene nanoribbons makes per- 
fect experiments and theory comparisons difficult. Both 
free-standing and supported ribbons are prone to 
distortions,^ ^ that affect also electronic structures; uni- 
axial stretching, for example, induces systematic changes 
in the fundamental energy gapsP^. Supported ribbons 
are easier to hold even,^but at the same time — because of 
support interactions — can have even more abrupt distor- 
tions such as loops and folds. ^. Ribbons can be affected 
by the support interactions also directly.^ 

Graphene nanoribbons have high claims for usage as 
sensors, spintronic devices, and ballistic transistors, but 
these claims originate often from simulations without de- 
fects, disorder or distort ions. Sure enough, distortions 
themselves, however complex, can be investigated by 
classical potentials and long, finite ribbons.^ ^ Bu t for 
electronic structure one needs a quantum approach ,1^^!^ 
which for finite ribbons would soon become expen- 
sive. Efficient twisting simulations hence require periodic 
boundary conditions — but beyond the standard transla- 
tional symmetry. 

In this Letter I show simulations of twisted zigzag 
(ZGNR) and armchair (AGNR) graphene nanoribbons, 
using density-functional tight-binding and revised peri- 
odic boundary conditions. This computationally efficient 
approach enabled simulations of infinitely long ribbons 
with continuously tunable twist. It turned out that, in 
terms of gap, density of states and optical spectrum, a 
twisted ribbon with a fixed length can be imitated by a 
fiat ribbon with a certain stretch. If the ribbon length 
is allowed to relax, however, the electronic structure re- 
mains almost unaffected by twisting; both of these phe- 
nomena can be explained in terms of average strain. 

The electronic structure modeling method was density- 
functional tight-binding, which is computationally ef- 
ficient and still captures the essentials of graphene's 
elastic and electronic properties. Spin was not in- 
cluded, leaving magnetic properties outside this work.^^ 
I adapted revised periodic boundary conditions to chiral 
symmetry, in the spirit of Ref. 15 , and used simulation 
cells as shown in Fig{l^; this approach is described in 

I simulated both bare-edge and hydrogen-passivated 



A^-AGNRs and A^-ZGNRs with A^ = 5 - 23, correspond- 
ing to widths W ^ 5 — 45 A. I characterize the twist by 
the dimensionless parameter rW, where r is the twist an- 
gle per unit length; ribbons with the same rW have the 
same strain at the edge and have the visual appearance 
of being equally twisted. Ribbons were twisted up to 
rW ~ 1, which corresponds to ribbons that are twisted 
full turns when ribbon lengths are 2ttW . Atom positions 
were optimized down to maximum force components of 
10~^ eV/A, and the simulation cell lengths were either 
fixed or fully relaxed .^^^^ 

Before turning attention to the electronic structure, it's 
illustrating to begin by inspecting geometry and energy. 
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FIG. 1: (Color online) a) Twisted 20- AGNR with rW = 0.25, 
together with an illustration of the simulation cell, b) Total 
energy per atom for twisted, hydrogen-passivated A-ZGNRs, 
either with fixed or relaxed unit cell length; only the narrowest 
ribbons deviate from the analytical estimates. Inset: energy 
can decrease upon twisting due to edge stress in unpassivated 
A-AGNRs, and induce spontaneous twisting. 
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The strain in the ribbon at a distance r from the axis is 

^(0=^. + ^, (1) 

where Ez is the axial strain (the unit cell strain). For 
6^=0 the edge strain is (rW)V8, or 12.5 % for rW = 1, 
and justifies the upper twist limit as a no-tearing con- 
dition. For a fixed length = 0) the elastic energy 
density ^ks'^ yields the total energy per atom the simple 
expression E = Ack/640{TW)^, where k = 25.5 eV/A^ 
is the in-plane modulus and Ac = 2.62 is the area 
per carbon atom; for a relaxed length the corresponding 
expression is E = Ack/lUOirW)"^. Figure [i]d, show- 
ing the energy for hydrogen-passivated ZGNRs, reveals 
that the strain really is the main component in energy, 
rendering these expressions fairly accurate despite their 
simplicity. The quartic energy dependence highlights the 
fiimsyness upon small twists {rW <C 1), suggesting usage 
in sensitive torsion balances — torsion constant is practi- 
cally zero. As soon as rW ~ 1, however, twisting rapidly 
becomes heavy. 

Hydrogen-passivated ribbons have no stress at the 
edge, but for ribbons with a compressive edge stress the 
situation is more involved. The inset in Fig. [TJd shows 
how twisting makes the energy for unpassivated AGNRs 
even to decrease, inferring spontaneous twisting. This 
phenomenon, which has been reported before,'^^ is sus- 
ceptible to edge passivation and reconstruction.f^^M^ In 
spontaneous twisting, however, energies are two orders of 
magnitude smaller than in forced twisting. 

Let us now turn attention to the results on electronic 
structure. Figure [2^ shows how the gaps of fixed-length, 
unpassivated A/'-AGNRs change upon twisting. Trends 
in the changes, just as in the gaps themselves (inset of 
Fig. [2|l), fall into three families that I define here as 

q = mod{N,3). (2) 

These trends agree with density-functional calculations of 
both finite and infinite ribbons.^ What is fascinating, 
though, is that for relaxed lengths these trends nearly 
vanish (Fig. [2]3). The three families are barely recogniz- 
able, substantial changes occurring only for the narrowest 
ribbons with W ^ 1 nm. The immediate question arises: 
how can one understand these differences? 

As it will turn out, the key for understanding these 
trends is Eq.Q. Fixed-length ribbons have £z = 0, 
and yield average strains of ^avg = (tW)^/24 across the 
ribbons. Now, Fig. compares three AGNRs that ei- 
ther are twisted with fixed length or are strained with 
£z = ^avg = (rW)^/24 and r = 0. The resulting sim- 
ilarity suggests the following interpretation: twisting by 
rW can be imitated through stretching by (rVK)^/24. The 
agreement is not perfect, and while it certainly could be 
improved by fitting, it would serve only little purpose, 
and merely give awkward dependencies on edge chirali- 
ties and passivations. 

Simple analytical expressions are also more powerful 
than dummy fitting parameters. Namely, combined with 
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FIG. 2: (Color online) a) Relative changes in energy gaps 
of unpassivated, fixed-length A'-AGNRs upon twisting. In- 
set: Absolute gaps for the corresponding flat ribbons; the q- 
families are defined in Eq.([2]). b) As above, but with relaxed 
length. Inset: Axial strains that result from length relax- 
ations follow the simple analytical expression. Line widths 
are proportional to ribbon widths. 

an axial strain e^, the average strain from Eq.Q is 

^Tavg = + (rVT) 2/24. (3) 

Minimizing the average strain l^avgl with respect to Sz 
for given rW minimizes also the total energy, inferring 
Ez = — (rW)^/24 and the trivial e^vg = 0. This provides 
a simple explanation to why energy gaps don't change in 
Fig. [2]3: with relaxed lengths the average strains across 
the ribbons are zero. 

The inset of Fig. shows that the relaxed length in- 
deed follows the analytical estimate accurately. With this 
estimate the position-dependent strain becomes e{r) = 
r^(r^/2 — iy^/24) — the strain is compressive at the axis, 
zero dit r ^ 0.29 • W, and (rW)^/12 at the edge. Length 
relaxation hence decreases the edge strain by two-thirds. 
The simulations agree surprisingly well with these simple 
expressions, especially in view of the noticeable anhar- 
monic effects due to > 10 % strains at the edge. 

The electronic structure is equivalent also beyond en- 
ergy gaps. Fig. [sJd compares the density of states and 
Fig.jSj) the imaginary part of the dielectric function, both 
for a 10-AGNR under four different conditions. In both 
figures a fiat ribbon and a relaxed- length twisted ribbon. 
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FIG. 3: (Color online) a) Relative changes in the energy gaps of three unpassivated A^-AGNRs, that were being either twisted 
with fixed length or only stretched, b) Density of states for a 10-AGNR that either was flat, was twisted to rW = 0.73 with 
relaxed length, was twisted to rW = 0.73 with fixed length, or was only stretched to Sz = (0.73)V24 = 0.022 (with r = 0). b) 
Imaginary part of the dielectric function under the same conditions as in b-panel. The dielectric function was calculated using 
random-phase approximation and polarization parallel to ribbon axis. 



as well as a fixed-length twisted ribbon and a stretched 
flat ribbon share the same features among themselves; 
the trend-determining factor is the average strain. 

Since most of the strain is at the edge, it is justified 
to ask whether twisting will bring about new localized 
electronic states at the edges. Closer analysis of the 
wave functions reveals, however, that no such localiza- 
tion takes place — also the the famous ZGNR edge states 
remained as before.!^ The changes in electronic structure, 
induced by twisting, proved to be smooth. 

To conclude, I have shown that changes in the elec- 
tronic structure of graphene nanoribbons upon twisting 
are determined by average strain. Because of the sim- 
plicity of this cause — although not all the results were 
displayed — same reasoning applies irrespective of ribbon 



chirality or edge passivation. Similar stretching-induced 
effects were recently reported also in carbon nanotubes 
under pure bending. Since flat, stretched ribbons have 
been investigated more, and since Eq. ([3| provides a 
direct mapping from stretched and twisted ribbons to 
stretched and flat ribbons, we thus have simple but pow- 
erful way to investigate ribbons under realistic experi- 
mental conditions. One should, however, be cautious 
not to take the twisting-stretching analogy too far — some 
properties, such as the ones related to magnetic flux,!^ 
may become fundamentally modified due to the uneven 
strain and the twisted geometry. 
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